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Abstract — This  work  examines  the  extension  of  a  recently 
developed  single  target  acoustic  localization  system  to  multiple 
targets.  Multiple  passive  sensors  distributed  in  a  region  of  space 
will  detect  the  acoustic  signal  of  multiple  gunfire  events  and 
measure  the  bearing  to  the  event  as  well  as  the  time  the  signal 
was  detected.  In  order  to  alleviate  the  data  association  problem, 
an  assignment-based  method  was  developed  in  order  to  associate 
measurements  at  a  fusion  center,  which  will  provide  estimates  of 
the  multiple  shooter  locations.  The  assignment  based  approach 
is  shown  to  localize  multiple  closely  spaced  shooters  with  a  high 
degree  of  accuracy. 

Keywords — acoustic  localization,  fusion,  passive  localization, 
assignment 


I.  Introduction 

This  work  extends  the  passive  gunfire-localization  system 
described  in  [5],  [6],  [14],  [15]  to  allow  for  localization  of 
multiple  shots  and/or  shooters. 

In  [14],  [15],  a  distinction  is  made  between  the  “native” 
measurements  (as  obtained  by  the  sensors)  and  the  “local 
estimates”  (based  on  the  native  measurements)  provided  by 
each  sensor  node  in  the  system.  The  native  measurements 
of  the  system  consist  of  bearing  (direction  of  arrival  — 
DOA)  measurements,  shockwave  DOA  measurements,  and 
time  difference  of  arrival  (TDOA)  between  the  times  of  each 
DOA  measurement.  The  bearing  measurements  are  provided 
by  measuring  the  DOA  of  the  muzzle  blast  of  a  gunshot,  while 
the  shockwave  DOA  is  a  measurement  of  the  shockwave  that 
emanates  from  a  supersonic  bullet  as  it  passes  the  sensor.  In 
the  original  system,  the  sensor  nodes  form  local  estimates  of 
the  range  and  the  direction  of  the  bullet  trajectory.  Accurate 
formation  of  the  local  estimates,  however,  depends  on  properly 
associating  the  bearing  and  shockwave  measurements  which 
pertain  to  the  same  gunfire  event.  Since  the  individual  sensor 
nodes  are  unlikely  to  correctly  associate  multiple  closely 
spaced  measurements,  the  remainder  of  this  work  will  not 
consider  local  estimates. 

Prior  research  in  acoustic  localization  of  small-arms  fire 
can  be  found  in  [1],  [3],  [5]-[13],  [17]-[19],  Some  of  these 
papers  use  “local  estimates”,  while  others  use  either  “native 
measurements”  or  only  time-of-arrival  (TOA)  or  TDOA  mea¬ 
surements.  Some  of  the  papers  [9],  [10],  [12],  [13]  employ 
a  bullet  deceleration  model,  while  others  assume  a  constant 
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velocity  bullet.  A  bullet  velocity  model  is  not  necessary  in 
this  work,  however,  since  shockwave  measurements  are  not 
used. 

In  this  work,  only  bearing  measurements  and  the  time  at 
which  the  bearing  measurement  is  detected  will  be  utilized 
for  multiple  shooter  localization.  The  framework  developed, 
however,  will  allow  inclusion  of  the  shockwave  measurements 
in  a  fairly  straightforward  manner  at  a  later  date.  The  fusion 
center  will  perform  the  association  via  the  generalized  mul¬ 
tidimensional  assignment  algorithm  (called  S-D  assignment 
in  [4]).  It  is  important  to  note  that  the  fusion  center  has 
no  knowledge  of  the  true  number  of  events.  The  assignment 
algorithm  will  look  for  the  overall  best  assignments  from  the 
available  measurements  at  each  sensor.  The  final  assignments 
will  then  undergo  a  threshold  test  on  the  number  of  non¬ 
dummy  (real)  measurements  included  in  each  assignment.  If 
the  number  of  non-dummy  measurements  in  an  assignment 
exceeds  the  threshold,  the  assignment  will  be  declared  as 
belonging  to  a  real  event.  The  remaining  events  will  be 
considered  false  events  by  the  fusion  center. 

Section  II  provides  the  details  of  the  problem  formulation 
and  introduces  the  notation  used.  Section  III  provides  simu¬ 
lation  results  for  a  specific  multiple  shooter  scenario.  Finally, 
Section  IV  will  provide  concluding  remarks. 


II.  Problem  Formulation 


Assume  there  are  Nt  (unknown  to  the  fusion  center) 
shooting  events  tpi,  i  =  1, . . . ,  NT  where 

V’i  =  [x*  tf!  (i) 

where  xt  =  [a is  the  location  (in  two-dimensional  space) 
of  the  z'th  shooting  event,  and  if  (“emission  time”)  is  the  time 
that  the  ?'th  shooting  event  occurred.  A  special  case  of  tjj0  will 
denote  a  false  event  (false  alarm). 


It  is  assumed  there  are  Ng  acoustic  sensors,  at  known 
locations 


(S1,...,SNs) 


SxN„  \ 

kJJ 


(2) 


observing  the  muzzle  blasts  of  the  shooting  events.  It  will  be 
assumed  that  each  sensor  receives  at  most  one  measurement 
for  each  true  gunfire  event  (i.e.,  there  is  assumed  to  be  no 
multipath  for  the  acoustic  signal,  however,  the  sensor  may  fail 
to  detect  the  muzzle  blast). 
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The  measurement  of  event  Vi  from  the  jth  sensor  is 

z3(i/>i)  =  h(ipi,  Sj)  +  Wj  (3) 

where 


'  tzi  L _ 

)2  +  (xi  -  SXj)2 J 

(4) 

where  0t  is  the  azimuth  (bearing)  from  sensor  j  to  the  location 
of  event  Vi>  t*  (“sensor  time”)  is  the  time  that  the  acoustic 
signal  from  the  muzzle  blast  of  event  Vi  reaches  sensor  j , 
and  c  is  the  propagation  speed  of  the  signal  (speed  of  sound 
through  air).  The  measurement  noise  of  sensor  j  is 


h(tpi,Sj) 


tan  y 

fi  +  W(y*  -  svi 


Wj  ~  TV  (0,  Rj) 


(5) 


where 


Rj  — 


’8,3 

0 


0 

7h 


(6) 


and  a |  •  and  <r2  j  are  the  variances  of  the  azimuth  and  time 
measurements,  respectively. 


Sensor  j  provides  a  set  of  Nj  measurements 

Zj  =  {z{,...,  4.},  j  =  l,...,Ns  (7) 

The  pdf  of  a  measurement  from  the  yth  sensor  is 

p(zJ  |  Vo)  =  p(0jli/>o)p(tjli>o)  (8) 

P(z3\ipi)  =  1 27t Rj | ~a/2  exp |  ^  [~  -  h(ipi,Sj)]'  RJ1 

■  [z3  -  h(ipi,Sj)]^,  i  =  Nt  (9) 

where  (8)  is  the  pdf  of  a  measurement  from  a  false  event 
( i  =  0,  i.e.,  a  false  alarm),  and  (9)  is  the  pdf  of  a  measurement 
from  a  true  event  ( i  0). 


The  lists  of  measurements  from  each  sensor  (7)  —  there 
are  Ns  lists,  i.e.,  the  dimension  of  the  assignment  is  Ns  — 
will  be  augmented  by  a  dummy  measurement  (indexed  by  0,  to 
represent  missed  detections).  The  TVg-tuple  of  measurements 
(representing  associated  measurements)  consisting  of  one  mea¬ 
surement  from  each  augmented  list  will  be  denoted 

ZWns  ={zl,zl,...,zl*s}  (10) 

where  £j  =  0, 1, . . . ,  Nj  represents  the  index  of  the  measure¬ 
ment  from  the  (augmented)  set  Z:l  which  is  included  in  the 
association. 1 


The  cost  of  a  given  association  ce1e2...eNs  will  be  given  by 
the  generalized  negative  log-likelihood  ratio 


C<1<2 


=  -  In 


A  (Ztit 


.IV  ) 


■  •^jvs  IVo) 


(ID 


where  VML  is  the  maximum  likelihood  (ML)  estimate  of 
the  shooting  event  pertaining  to  the  particular  association  of 
measurements  Z(1c2,,jNs. 


'Recall  that  tj  =  0  represents  the  dummy  measurement,  so  (10)  need  not 
contain  Ng  “real”  measurements,  i.e.,  missed  detections  are  allowed  in  the 
association. 


Assuming  the  measurements  are  independent  across  sen¬ 
sors,  the  likelihood  function  of  a  particular  shooting  event  Vi> 
for  the  association  is 

Ns 
3= 1 

/  \  “(D 

•  [PdjP(z3e.  |  Vi)J  i  =  l,...,NT 

(12) 

where  /-)/,  is  the  probability  of  detection  for  sensor  j,  and  the 
indicator  function  u{£j)  is 


(  ,a  r°  if v-  =  o 

[  j>  ~  1  1  otherwise 


03) 


The  likelihood  function  of  the  false  shooting  event  Vo  is 

-r-T  /  \u(tj) 

A {Ztle2...eNs  |Vo)  =  (p(^|Vo)p(i|JVo)J  (14) 

3=1 

The  pdf  of  the  muzzle  blast  DOA  measurement  in  the  case  of 
a  false  shooting  event  is  assumed  to  be 

P(0\ipo)  =  (15) 

i.e.,  the  false  bearing  measurements  are  assumed  to  be  uni¬ 
formly  distributed.  The  pdf  of  the  detection  time  of  a  false 
shooting  event  is 

P(^lVo)  =  ^7  (16) 

where  W  is  the  time  window  for  the  events  under  consid¬ 
eration.  For  events  and  measurements  which  are  separated 
significantly  in  time,  there  is  no  data  association  ambiguity, 
so  it  is  assumed  that  only  measurements  falling  within  the 
time  windows  W  need  to  be  associated. 


The  final  cost  of  association  is  now 

Chi 2-Lvs  =  In  A  {zeit2...lNs  |Vo)  -  In  A  {Zt^2...tNs  |Vi) 

Ns  Ns 

=  -5>(^)ln(27 rlVPd.)  - 

3=1  3=1 

.  NS 

+  -Y,  uVi)  ln(|27Ti?j|)  [z3  -  h(tpi,  Si)] '  RJ1 

3=1 

■  [z3  -  Sj)\  (17) 


^(l-«(^))ln(l-Pdj) 


A.  Calculation  of  VMi  and  Final  Assignments 

Since  the  cost  of  a  candidate  association  (11)  requires 
the  ML  estimate  of  the  shooting  event  pertaining  to  that 
association,  the  Iterated  Least  Squares  (ILS)  estimator  [2] 
will  be  used  to  provide  the  estimate.  This  is  identical  to  the 
estimation  approach  used  in  [14],  [15],  which  was  shown 
to  be  statistically  efficient.  The  ML  estimate(s)  associated 
with  the  final  assignment(s)  will  also  be  provided  as  the 
estimate(s)  of  the  target  location(s)  at  the  conclusion  of  the 
association  algorithm.  Then  the  fusion  center  will  output  a 
list  of  measurement  associations,  as  well  as  a  list  of  target 
locations  (and  estimates  of  the  time  the  shot  was  taken,  if  that 
is  desired). 
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TABLE  I:  Fraction  of  Pure  and  Complete  (PC),  Pure  and 
Incomplete  (PI),  Impure  2  (12),  and  Impure  1  (II)  assignments 
(in  1000  Monte  Carlo  runs). 


•  Sensor 

★  Target 


* 


S4 

:  • 


S5 

:  • 


S2 


S3 

■■■•■ 


-50  -25  0  25  50  75  100  125  150  175  200  225  250 

Fig.  1:  Overhead  view  of  the  multiple  shooter  scenario. 


Pd 

Nfa 

r 

PC 

PI 

12 

11 

0.90 

1.28 

2 

0.791 

0.179 

0.000 

0.030 

3 

0.909 

0.090 

0.001 

0.000 

0.95 

1.28 

2 

0.934 

0.039 

0.000 

0.027 

3 

0.974 

0.026 

0.000 

0.000 

0.99 

1.26 

2 

0.990 

0.000 

0.000 

0.010 

3 

1.000 

0.001 

0.000 

0.000 

0.90 

2.49 

2 

0.728 

0.185 

0.001 

0.086 

3 

0.890 

0.109 

0.001 

0.000 

0.95 

2.50 

2 

0.842 

0.071 

0.000 

0.087 

3 

0.945 

0.055 

0.001 

0.000 

0.99 

2.50 

2 

0.954 

0.002 

0.000 

0.044 

3 

0.998 

0.002 

0.000 

0.000 

The  necessary  terms  for  the  Jacobian  of  (4), 
perform  the  ILS  algorithm,  are 

90 j  yi  Syj 

dxi  (xi  -  Sx. )2  +  ( yi  -  SVj)2 

dO j  Xi  dx,j 

dyt  { Xi  -  SXj )2  +  {yi  -  Sy.)2 


dtj  _ _ Xi  -  SXj _ 

dxi  c^/{xi  -  S x j ) 2  +  {yi  -  Syj)2 

= _ Vi  -  Syj _ 

dyi  CyJ{xi  -  SXj)2  +  {yi  -  Sy.)2 
dtd 

-—d-  =  1 


in  order  to 

(18) 

(19) 

(20) 
(21) 
(22) 
(23) 


Each  sensor  has  an  identical  probability  of  detection  /)/,  as 
well  as  an  identical  spatial  density  of  false  alarms.  The  false 
alarms  are  assumed  distributed  uniformly  in  [0,  27t]  for  the 
DOA  measurements  and  uniformly  in  [0,  W]  for  the  detection 
time  (the  time  window  W  in  this  case  is  set  to  1  s). 

Once  the  assignment  algorithm  is  finished,  any  assignments 
with  fewer  than  r  non-dummy  measurements  are  discarded 
(since  a  single  non-dummy  measurement  cannot  provide  a 
valid  gunfire  event  estimate,  r  will  generally  be  either  2  or  3). 
The  assignments  passing  this  threshold  test  are  then  analyzed 
as  follows,  for  varying  levels  of  Pd,  r,  and  Nfa  (the  average 
number  of  false  alarms  from  each  sensor  in  the  time  window 
W  =  1  s). 

Table  I  shows  the  fraction  of  the  total  assignments  that  fall 
into  four  categories 


Once  the  costs  (11)  of  the  candidate  associations  are  calcu¬ 
lated,  the  generalized  multidimensional  assignment  algorithm 
of  [4]  is  performed  to  form  the  final  measurement  associations. 


(PC)  Pure  and  Complete:  the  non-dummy  measurements  of 
the  assignment  all  come  from  the  same,  real  (i.e.,  not 
a  false  alarm)  source,  and  all  the  measurements  of  the 
event  are  included. 


III.  Simulation  Results 

In  the  following  simulations  there  are  two  true  gunfire 
events 


(PI)  Pure  and  Incomplete:  the  non-dummy  measurements 
all  come  from  the  same,  real  source,  but  not  all  the 
measurements  of  the  event  are  included. 


=  [0  200  0.2]' ,  =  [2  202  0.25]'  (24) 


separated  by  2  m  (in  each  spatial  dimension)  and  50  ms. 


There  are  five  sensors  located  at 


Si 


127 
107j  ’ 

S4  = 


S2  = 

136 
68  ’ 


20 

22j  5 

S5  = 


S3  = 

182 

59 


90 
0  ’ 


(25) 


The  overhead  view  of  this  scenario  is  shown  in  Figure  1. 
The  standard  deviations  of  the  measurement  noises  are 


ogj  =  1  mrad  j  =  1, . . . ,  5  (26) 

<Jtj  =  lms  j  =  1,  ...,5  (27) 


(12)  Impure  2:  at  least  two  non-dummy  measurements  from 
the  same,  real  source  (still  useful). 

(II)  Impure  1:  no  two  non-dummy  measurements  are  from 
the  same,  real  source  (ghost). 

When  comparing  the  results  in  Table  I  for  r  =  2  and  r  =  3, 
note  that  the  number  of  total  assignments  will  be  higher  for 
the  lower  threshold.  What  this  means  is  that  even  though  the 
fraction  of  PC  assignments  is  lower  for  r  =  2,  the  number 
of  assignments  which  are  PC  may  be  similar.  What  Table 
I  demonstrates  is  the  effect  of  the  threshold  on  the  impure 
assignments,  as  well  as  the  effect  that  Pd  has  on  the  purity 
of  the  final  assignments.  Obviously,  it  is  desirable  to  have  the 
purest  possible  assignments  to  ensure  the  final  estimate  is  as 
accurate  as  possible. 
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TABLE  II:  Probability  of  detecting  true  events,  probability  of 
split  assignments,  and  average  number  of  false  events  detected 
(in  1000  Monte  Carlo  runs). 


Pd 

Nfa 

T 

Event 

Pn, 

Ps, 

Nfd 

o 

1 

0.914 

0.086 

0.067 

0.90 

1.28 

2, 

2 

0.918 

0.081 

3 

1 

0.957 

0.001 

0.000 

2 

0.959 

0.000 

2 

1 

0.990 

0.010 

0.056 

0.95 

1.28 

2 

0.988 

0.012 

q 

1 

0.995 

0.000 

0.000 

o 

2 

0.994 

0.000 

1 

1.000 

0.000 

0.020 

0.99 

1.26 

A 

2 

1.000 

0.000 

q 

1 

1.000 

0.000 

0.000 

D 

2 

1.000 

0.000 

1 

0.919 

0.081 

0.204 

0.90 

2.49 

z 

2 

0.917 

0.083 

3 

1 

0.956 

0.000 

0.000 

2 

0.961 

0.002 

2 

1 

0.982 

0.018 

0.193 

0.95 

2.50 

2 

0.983 

0.017 

'i 

1 

0.988 

0.001 

0.000 

o 

2 

0.991 

0.000 

1 

1.000 

0.000 
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Table  II  shows  a  summary  of  the  event  detection  results 
of  the  assignments.  The  probability  of  detecting  true  event 
ipi,  Pnt ,  is  calculated  as  the  fraction  of  runs  where  only  a 
single  assignment  contained  two  or  more  measurements  of  the 
true  event  ipi.  The  probability  of  having  the  measurements 
of  event  ipi  split  over  more  than  one  assignment  is  denoted 
PSi  (to  counteract  the  effect  of  split  assignments,  a  test  could 
be  devised  for  merging  assignments).  The  probability  of  not 
detecting  event  ipi  (miss)  could  then  be  calculated  as 

PMi=l-(PDi+Psi)  (28) 

Finally,  the  average  number  of  “false  events”  detected  per  run 
is  denoted  NFd-  An  assignment  which  gives  rise  to  a  “false 
event”  is  defined  as  one  in  which  one  or  more  false  alarms  are 
associated  to  no  more  than  one  real  measurement  from  any 
actual  gunfire  event. 

In  this  case,  the  lower  threshold  can  be  seen  to  allow 
detections  of  false  events,  while  the  higher  threshold  (for 
this  scenario)  completely  eliminates  detections  of  false  events. 
The  lower  threshold,  however,  provides  a  nearly  zero  miss 
probability  Pm,  ,  while  the  higher  threshold  will  miss  the 
detection  of  the  true  events  approximately  4%  of  the  time  for 
Pd  =  0.9.  For  higher  individual  sensor  detection  probabilities, 
however,  the  miss  probability  is  quite  low,  even  for  the  higher 
threshold. 

Figures  2a-7b  show  the  final  estimation  results  of  the 


various  configurations  shown  in  Tables  I  and  II.  The  figures 
with  t  =  3  show  the  high  accuracy  of  the  final  associations, 
with  distinct  groupings  around  each  of  the  two  true  events.  For 
the  figures  with  r  =  2,  the  detection  of  false  events  results 
in  estimates  which  are  far  from  the  true  events,  making  the 
accuracy  of  the  true  event  estimates  difficult  to  judge  at  that 
scale.  For  this  reason,  inset  plots  are  included  with  a  scale 
that  matches  the  corresponding  r  =  3  figure.  These  inset  plots 
show  similar  estimate  results,  however,  for  the  lower  Pd  cases 
there  is  a  larger  spread  (due  to  split  assignments  and/or  false 
events). 

The  average  computation  time  for  the  generalized  multidi¬ 
mensional  assignment  algorithm  to  complete  one  of  the  1000 
Monte  Carlo  runs  is  0.35  s.  This  computation  was  performed 
on  a  2.66GHz  Intel  Core2  Windows  machine,  running  Mat- 
lab  2013b,  with  the  generalized  multidimensional  assignment 
running  as  a  “mex”  file  written  in  C++. 

IV.  Conclusions 

This  work  extends  existing  methods  of  gunfire  localization 
systems  which  are  meant  to  localize  a  single  gunfire  event,  so 
that  the  localization  system  is  capable  of  localizing  multiple 
events  (i.e.,  multiple  shooters  and/or  multiple  shots).  This 
particular  system  utilizes  only  bearing  and  time  measurements, 
and  has  no  knowledge  of  the  true  number  of  events.  By 
leveraging  existing  assignment  methods  to  solve  the  data  asso¬ 
ciation  problem  in  this  multiple  shooter  case,  the  assignment 
algorithm  is  able  to  very  accurately  localize  multiple  closely 
spaced  gunfire  events.  A  simple  threshold  test  on  the  number 
of  non-dummy  measurements  in  a  final  assignment  serves  to 
eliminate  assignments  that  would  lead  to  false  events.  A  lower 
threshold  will  give  rise  to  false  event  detections  as  well  as 
“split  assignments”  (which  will  necessitate  further  tests  for 
merging  assignments).  A  higher  threshold  will  result  in  a 
slightly  higher  missed  detection  rate  of  true  events,  but  will 
suppress  false  events  and  split  assignments.  Future  work  will 
extend  this  method  to  incorporate  shockwave  measurements, 
examine  the  impact  of  possible  multipath  propagation  of  the 
acoustic  signals,  and  examine  different  methods  of  sequential 
assignment  as  alternatives  to  full  multidimensional  assignment. 
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